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Abstract 


We carried out the first multi-wavelength (optical/UV and X-ray) photometric reverberation mapping of a tidal 
disruption flare (TDF) ASASSN-14li. We find that its X-ray variations are correlated with and lag the optical/UV 
fluctuations by 32 + 4 days. Based on the direction and the magnitude of the X-ray time lag, we rule out X-ray 
reprocessing and direct emission from a standard circular thin disk as the dominant source of its optical/UV 
emission. The lag magnitude also rules out an AGN disk-driven instability as the origin of ASASSN-14li and thus 
strongly supports the tidal disruption picture for this event and similar objects. We suggest that the majority of the 
optical/UV emission likely originates from debris stream self-interactions. Perturbations at the self-interaction sites 
produce optical/UV variability and travel down to the black hole where they modulate the X-rays. The time lag 
between the optical/UV and the X-rays variations thus correspond to the time taken by these fluctuations to travel 
from the self-interaction site to close to the black hole. We further discuss these time lags within the context of the 
three variants of the self-interaction model. High-cadence monitoring observations of future TDFs will be sensitive 
enough to detect these echoes and would allow us to establish the origin of optical/UV emission in TDFs in 
general. 
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1. Introduction 


When a star comes sufficiently close to a massive black hole 
(1075 Mo) such that tidal forces exceed its self-gravity, it will 
be torn apart to cause a stellar tidal disruption flare (TDF; 
Hills 1975; Rees 1988). These events not only have the 
potential to uncover a large population of hidden massive black 
holes but could also enable us to understand how jets are 
launched and disks are formed. 

One of the recent controversies about TDFs is, where does 
the majority of the optical/UV emission originate from in 
TDFs? (Piran et al. 2015 (P15); Roth et al. 2016; Dai 
et al. 2015; etc.). In some of the best-studied events, the 
optical/UV blackbody radius is much larger than the tidal 
radius, and the temperature remains roughly constant over 
much of the evolution (e.g., Chornock et al. 2014; Holoien 
et al. 2016b). Because TDFs are basically accreting super- 
massive black holes, the natural inclination is to expect the 
same mechanisms operating in AGN, namely, X-ray reproces- 
sing, thermal emission from the inner radii of the accretion 
disk, etc., to also operate in TDFs. Here, we use publicly 
available multi-wavelength (optical/UV and X-ray) data from a 
recent TDF, ASASSN-14li, to constrain the origin of its 
optical/UV emission. 

ASASSN-14li was discovered as an optical transient by the 
All-Sky Automated Survey for SuperNovae (ASASSN; 
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Shappee et al. 2014) on 2014 November 11 (Holoien 
et al. 2016a, hereafter H16). Based on its spatial coincidence 
with its host galaxy’s nucleus, its peak luminosity of roughly 
10“ergs~', and its blue optical spectrum with broad and 
transient Ho and He emission lines, it has been categorized as 
an event caused by the tidal disruption of a star by a 
supermassive black hole (H16). We carried out a photometric 
reverberation analysis of ASASSN-14li to find that the 
variations in its X-ray bandpass are correlated with and lag 
the optical/UV fluctuations by 32 days. We describe the data 
and our analysis in Sections 2 and 3 and discuss the 
implications of our result for the likely origin of the optical/ 
UV emission from this source in Section 4. 


2. Optical/UV and X-Ray Observations 


The multi-wavelength data used in this Letter were acquired 
from two different facilities: the Swift satellite (Gehrels 
et al. 2004) and the Las Cumbres Observatory Global 
Telescope (LCOGT) network (Brown et al. 2013). Swift's 
X-Ray Telescope (XRT; Burrows et al. 2005) provided the 
X-ray (0.3-10 keV) data, while its UV Optical Telescope 
(UVOT; Roming et al. 2005) provided the UV data in the 
UVW2, UVM2, and the UVWI filters with centroid wave- 
lengths of 1928, 2246, and 2600 A, respectively (Poole 
et al. 2008). The UVOT also facilitated optical data in the U, 
B, and the V bands with centroid wavelengths of 3465, 4392, 
and 5468 A, respectively (Poole et al. 2008). However, because 
the B- and the V-band optical observations were poorly 
sampled for the second half of Swiffs monitoring campaign 
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(see Figure 2 of H16), we also used published LCOGT's 
g-band optical data (centroid wavelength of 4770 A). We re- 
analyzed all of the X-ray and the UVOT data from Swift. Our 
reduction and analysis procedures are described in detail below. 


2.1. X-Ray Data Analysis 


Swift started monitoring the TDF ASASSN-14li roughly 
eight days after its discovery on MJD 56983.6 (H16). It 
observed the source for 1-3 ks roughly once every three days 
for the first 2:260 days of the campaign. However, after this 
time the target became Sun-constrained to Swift and the 
monitoring campaign suffered from longer data gaps. There- 
fore, we limited our analysis to the first 270 days of post- 
outburst data at all wavelengths. 

After extracting the clean eventlists from the raw Level-1 
XRT products, we accounted for pile-up (Miller et al. 2015) by 
extracting X-ray (0.3-10.0 keV) source events from an annular 
region (rather than a circular region) centered on the source's 
centroid. As recommended by the XRT data analysis guide,'! 
we modeled the PSF in each observation to estimate an inner 
exclusion radius. The inner exclusion radius of the source 
extraction region varied between 0" and 15" while the outer 
radius was fixed at 50". The background count rates were 
extracted from annuli centered on the source with inner and 
outer radii of 70" and 250", respectively. Thus, we extracted a 
pile-up and background-corrected source count rate from each 
Swift observation. We also corrected for bad pixels followin 
the procedure outlined by the Swift data analysis guide." 
Because ASASSN-14li’s X-ray energy spectrum remains 
roughly constant during the first 250 days of the outburst (see 
Figure 3 of Miller et al. 2015 and Table S5 of Brown et al. 
2017), the count rate serves as a good indicator of its 
intrinsic flux. 


2.2. UVOT Data Analysis 


The UVOT (Roming et al. 2005) on board Swift began 
observing ASASSN-14li at 10:08 UT on 2014 November 30. 
We downloaded the data from the HEASARC archive." In 
each of the six broadband filters (UVW2, UWM2, UVWI1, U, 
B, and V), we corrected the image astrometry and stacked 
frames on a per-visit basis. 

We performed photometry on the stacked images using a 3” 
aperture and corrected to the standard UVOT photometric 
system (Poole et al. 2008; Breeveld et al. 2011). To account for 
contamination from the host underlying a host galaxy, a 
particularly important contribution to the observed emission in 
the redder filters, we utilize the coincidence-loss corrected host 
flux estimates from Miller et al. (2015; see also H16 and van 
Velzen et al. 2016). All the optical/UV light curves are 
consistent with H16. The cross-correlation results below are 
independent of host-subtraction. 


3. Cross-correlation Analysis 


In order to search for and quantify any correlation between 
the X-ray and the optical/UV variations, we employed the 
interpolated cross-correlation function (ICCF) methodology as 
described by Peterson et al. (2004). We first de-trended the 


H http:/ /www.swift.ac.uk/analysis/xrt/pileup.php 
'2 http: //www.swift.ac.uk /analysis/xrt/Iccorr.php 
13 See http: / /heasarc.gsfc.nasa.gov /cgi-bin/W3Browse/swift.pl. 
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X-ray light curve with a bending power-law model (red curve 
in the top left panel of Figure 1). Each of the optical/UV light 
curves were then de-trended with a power-law decay model. A 
sample fit for the UVW1 data is shown as a dashed red curve in 
the top middle panel of Figure 1. The X-ray residual was then 
cross-correlated with each of the seven optical/UV residuals. A 
sample ICCF between the X-ray and the UVWI data is shown 
in the top right panel of Figure 1, while the rest of the ICCFs 
are shown in Figures 2 and 3. It is evident that the CCFs peak 
between 25 and 35 days indicating that the X-ray variations lag 
the optical/UV fluctuations by 25-35 days. 

In order to confirm these CCFs, we also visually inspected 
them and compared the variability features in the X-ray data 
with the features in the optical/UV residual light curves. These 
are shown in the bottom right panel of Figure 1 and the right 
panels of Figures 2 and 3. In each of these panels, common 
features in both the light curves are evident. Some of these 
CCFs also show evidence for an anti-correlation at zero lag. 
However, we do not consider it to be real as visual inspection 
of the light curve residuals did not reveal a strong evidence for 
an anti-correlation. 

The X-ray and the lead-corrected optical/UV residual light 
curves in Figures 1—3 suggest that the correlation functions are 
driven by three X-ray variability features: two broad peaks (one 
around day 50 and the second around day 110) and one trough 
around day 80 in the X-ray data (see the top right panel of 
Figure 2). Correspondingly, these features are also present in 
the optical/UV residual light curves but offset by ~—32 days 
(see the top right panel of Figure 2). 

In order to estimate the statistical significance of the 
correlations, we estimated the chance probability that the 
X-ray light curve would produce these same variability features 
evident in a given optical/UV light curve. For this purpose, we 
simulated 10? white noise light curves to mimic the X-ray 
residual light curve. We constructed these synthetic X-ray light 
curves by drawing from a Gaussian distribution with a mean of 
zero and a standard deviation equal to the standard deviation of 
the observed X-ray residual light curve. We then evaluated 
simulated ICCFs by cross-correlating each of these simulated 
X-ray residual light curves with each of the optical/UV 
residual light curves. This way—for each X-ray—optical/UV 
ICCF—we built a distribution of the 10° correlation values at 
each lag. Using these distributions, we extracted the 9596 (blue) 
and the 99% (magenta) confidence contours shown in the top 
right panel of Figure 1 (see also Figures 2 and 3). We then 
estimated the uncertainty in the peak and the centroid of the 
ICCFs using the Random Subset Selection and the Flux 
Randomization procedures as described by Peterson et al. 
(1998, 2004). The resulting cross-correlation centroid and peak 
distributions for the case of X-ray versus UVW1 ICCF are 
shown in the bottom left and the bottom middle panels of 
Figure 1, respectively. 

Using the binomial distribution formula, we also estimated 
the global statistical significance of the X-ray variations lagging 
the optical/UV fluctuations to be >4.4c. This was estimated 
based on the fact that five out of seven CCFs have the 
correlation peaks detected at greater than the 9596 confidence 
(see, for example, Pasham et al. 2015 and Tombesi et al. 2010 
within the context of timing and spectroscopy, respectively). 

We also evaluated all the X-ray—optical/UV CCFs using the 
discrete cross-correlation function (DCF) algorithm as 
described by Edelson & Krolik (1988). We find that the 


Pasham et al. 


THE ASTROPHYSICAL JOURNAL LETTERS, 837:L30 (9pp), 2017 March 10 


X-ray Light Curve UVW1 Light Curve 


0.6 Af 
m 0.5 gH V 
2 E Z 
8 03 B2 5 
V £ S 
© 
É 0.2 p^ E 
[za 31 QU 
0.1 E 
0.0 : 
50 100 150 200 250 50 100 150 200 250 -20 0 20 40 
Time (in days) Time (in days) X-ray Lag (in days) 
CCF Centroid Distribution CCF Peak Distribution X-ray Vs UVW1 
3000 ; 
8. 2500 g 3000 AN 
a S 5 
$ 2000 AN 8 
E 3 2000 EB 
m H 
* 1500 $ E: 
B 9 1500 5 
vo I^] iss) 
È 1000 2 1000 x 
3 E = 
Z 500 Z 500 rs 
0 0 
10 20 30 40 50 60 10 20 30 40 50 60 50 100 150 200 250 


Lag Centroid (in days) Time (in days) 


Peak Lag (in days) 


Figure 1. Top left: ASASSN-14li’s observed X-ray light curve (black data points) along with the best-fit bending power-law model (dashed red). Top middle: 
ASASSN-14li's observed UVWI light curve (blue) along with the best-fit decaying power-law model (dashed red). All the best-fit model parameters can be found in 
the data behind this figure. Top right: interpolated cross-correlation function (ICCF) between the de-trended X-ray and the UVWI light curves. The 95% (blue) and the 
99% (magenta) white noise statistical confidence contours are also shown. The solid vertical red line is the median of the CCF’s centroid distribution (3572 days) 
shown in the bottom left panel. The lo error bars are shown as the dashed vertical red lines (same as in the bottom left panel). Bottom middle: ICCF’s peak 
distribution. The solid vertical line is the median while the dashed lines are the lo error bars. Bottom right: we compare the variability features in the X-ray and the 
UVW! light curves. ASASSN-14li's UVWI light curve is offset by 35 days (blue). The black data points show the X-ray light curve interpolated at the lead-corrected 
UVW1 epochs. Both the light curves were de-trended by subtracting a smooth function (top panels), leaving only the variability features. The solid curves are a 
running average of five neighboring points, to guide the eye. All the X-ray and optical/UV photometric measurements can be found in the data behind this figure. 


(The data used to create this figure are available.) 


resulting DCFs are consistent with the ICCFs shown in 
Figures 1-3. 

Because a CCF is simply a convolution of the autocorrela- 
tion function (ACF) with a transfer function, a lag between the 
X-ray and the optical/UV light curve is real only if it arises 
from the transfer function and not from the X-ray or the 
optical/UV ACF. To rule out an ACF origin for the lag we 
extracted the X-ray and the optical/UV ACFs to find that they 
do not have any statistically significant lag features except at 
zero lag. 


4. Discussion 


First, we demonstrate that the measured optical/UV—X-ray 
lag in ASASSN-14li implies that none of the three primary 
mechanisms ^ that can drive X-ray variations to correlate with 
the optical/UV changes in AGN (e.g., Edelson et al. 2015; 
McHardy et al. 2016) are at play here. 


14 Other mechanisms similar to stellar-mass black hole binaries may also play 
a role (e.g., Gandhi et al. 2010 and references therein), but on timescales much 
shorter than observed here. 


A. In the first mechanism, X-rays from close to the black 
hole scatter off material in the outer regions of the 
accretion disk or any other surrounding medium, lose 
energy, and get reprocessed into lower-energy optical/ 
UV photons (Edelson et al. 2015). In this case, the X-ray 
variations would lead the optical changes by a few 
days to a few tens of days depending on the distance 
to the outer disk (Morgan et al. 2010)/reprocessing 
medium. 

B. However, if a significant fraction of the optical/UV 
emission originates directly from the accretion disk as a 
thermal blackbody, then two kinds of correlations are 
possible: 

(1) If the majority of the X-rays are produced in a corona 
that is powered by Compton up-scattering of near-UV 
seed photons from the thermal inner disk (Reynolds & 
Nowak 2003), then the X-ray variations would lag the 
near-UV fluctuations by a fraction of a day to a few 
days depending on the size of the corona and the inner 
disk (which are determined by the black hole mass 
(Arévalo et al. 2005) and accretion rate). 
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Figure 2. Left panels: X-ray vs. optical/UV cross-correlation functions (CCFs). The two cross-correlated light curves are indicated at the top of each panel. The blue 
and the magenta curves are the 95% and the 99% confidence contours, respectively. The solid vertical red line is the median of the CCF’s centroid distribution, while 
the dashed vertical lines are the 1c error bars on the median. These values for the UVW2, UVM2, and the UVW1 ICCFs are 26122, 3075, and 3279 days, respectively. 
Right panels: we compare the variability features in the X-ray light curve (black) with the optical/UV data (blue). The optical/UV light curves were offset by their 
corresponding lead times. The solid curves are the running average of five neighboring points to guide the eye. Common features in both bands are evident in all the 


cases. 
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Figure 3. Same as Figure 2, but the median of the CCF's centroid distribution for the SDSS-g, B, and the V band ICCFs are 371%, 3217o, and 31718 days, 
respectively. Top right: because the SDSS-g band light curve was better sampled than the X-ray data, we offset the X-ray light curve by the lag implied from the CCF 
and interpolated the SDSS-g onto the X-ray epochs. 
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(2) Also, in this scenario, optical/UV fluctuations due to 
accretion rate perturbations could propagate inward on 
a timescale corresponding to the local viscous time- 
scale (Arévalo et al. 2008). Then the optical/UV 
emission leads the X-rays by a few tens to millions of 
days (Breedt et al. 2009) again depending on the black 
hole mass and the accretion rate (Shakura & 
Sunyaev 1973; Figure 3). 


It is evident from Figures 1—3 that the optical/UV emission 
"leads" the X-rays by 32 + 4days. This rules out X-ray 
reprocessing as the dominant source of the variable optical/UV 
emission. Furthermore, assuming a black hole mass of 109? 0.6 
Mo as derived from its host galaxy’s bulge luminosity and the 
bulge's stellar velocity dispersion (H16, van Velzen 
et al. 2016), the seed photon scenario (case (1) above) is also 
unlikely as it would result in UV lead times of only a few 
thousands to a few tens of thousands of seconds. Below, we 
show that direct optical/UV emission from a standard thin disk 
can also be ruled out. 

ASASSN-l4li's observed peak bolometric luminosity of 
roughly 10^ erg s™' (H16) implies an accretion rate of 0.2707 
Mgaa (where Mgaa is the Eddington accretion rate). Such a sub- 
Eddington rate implies that if an accretion disk formed quickly 
after the disruption it can be described by a geometrically thin, 
optically thick disk model of Shakura & Sunyaev (1973). The 
observed lags between the optical, the UV, and the X-ray 
radiation, if coming from the viscous propagation of accretion rate 
fluctuations in such a thin disk, require, for possible black hole and 
disk parameters of ASASSN-14li, that the optical/UV emission 
come from within a radius of 50 R, (top left panel of Figure 4; R, 
is the gravitational radius). If emitted from outside this radius, the 
viscous timescales, and therefore the lags, would be much longer 
than the observed lags. On the other hand, the thin disk around a 
1095-056 M. black hole of the appropriate luminosity emits most 
of its optical and UV radiation from outside a radius of 500 R, (top 
right panel of Figure 4). As a result, the expected lags due to 
viscous propagation of perturbation in a thin disk are orders of 
magnitude longer than the observed ones (bottom panel of 
Figure 4). This discrepancy rules out the standard, thin, circular 
disk solution. More importantly, this argues strongly that 
ASASSN-14li is not an accretion disk-driven AGN flare. 

The slim disk accretion disk model (Abramowicz et al. 1988), 
describing super-Eddington accretion flows, predicts a thicker 
accretion flow that can produce shorter viscous times compar- 
able with ASASSN-14li’s lags. However, ASASSN-14li is sub- 
Eddington and in this limit the slim disk model reduces to the 
fiducial thin disk solution (Strubbe & Quataert 2009). Another 
possibility is that a circular disk of a different nature, thick and 
very hot (e.g., Coughlin & Begelman 2014), forms. Such a disk 
can form in principle even for sub-Eddington accretion rates. 
However, no general solution for such a mode exists, and 
therefore in this picture the lag magnitude cannot be used to 
constrain the flow parameters. Finally, a truncated disk model— 
where the optical/UV are produced by the Rayleigh-Jeans end 
of an X-ray multi-color blackbody emission—cannot explain the 
observed high optical/UV luminosities (H16; but see below). 

Many recent numerical studies of tidal disruption events (e.g., 
Guillochon et al. 2014, G14; Shiokawa et al. 2015; Hayasaki 
et al. 2016; Bonnerot et al. 2017, hereafter B17) have shown that 
the infalling stellar debris stream will undergo self-interactions 
because of relativistic apsidal precession (see, for example, 
Figure 2 of B17). While these studies differ in their predictions 
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for the evolution of the debris stream following the first self- 
interaction, they all agree that the role of self-interactions /shocks 
is to facilitate the process of circularization by removal of 
angular momentum. For instance, Shiokawa et al. (2015) and 
P15 argue that following the first stream self-interaction at the 
apocenter of the most bound debris, material spreads and rushes 
toward the black hole on the free-fall timescale. On the other 
hand, B17 suggested a discrete variant of the self-interaction 
model where the stream undergoes successive self-interactions— 
without much spreading in each interaction—and can eventually 
settle into a small (~ a few tens of gravitational radii) accretion 
disk. Alternately, the debris can instead lose only a small fraction 
of its energy in each self-interaction and follow more bound 
elliptical orbits. This would ultimately lead to the closing of the 
gap between the apocenter and the black hole, forming an 
elliptical accretion disk (G14). 

The lags of the observed magnitudes can originate in all the 
three models (see Figure 5 for schematics). In the shock-at- 
apocenter model, the collisions of streams are likely to result in 
perturbations of the, otherwise relatively uniform and thin, tidal 
stream of gas, and subsequently vary the fraction of the tidal 
debris getting close to the black hole and being trapped or 
directly accreted there (e.g., Sadowski & Narayan 2015). Such 
an inner small-scale accretion flow will be hot and will be 
modulated at the rate at which tidal debris—already affected by 
the self-interaction at large radii—Treturns to very close orbits. 
Therefore, heating up a clump of gas in the dissipation/ 
interaction region resulting in the optical/UV emission will be 
followed by modulation of the energetic X-ray radiation 
coming from a region close to the black hole. Furthermore, 
under such circumstances, X-rays variations are expected to lag 
the optical/UV fluctuations by the infall time from the stream— 
stream interaction region. This lag is roughly equal to half the 
orbital time of the debris orbit. For example, the shortest 
expected lag in this model would correspond to half the orbital 
period of the most bound orbit. This is roughly 11 days for a 
1095 Mo black hole (using Equation (4) of P15) and is in 
agreement with the lags identified in this work. It is also 
plausible that fluctuations propagate on a thermal timescale that 
would be slower than the free-fall timescale. 

Bonnerot et al. (2017) argue that in the absence of strong 
magnetic stresses, successive self-interactions can lead to the 
formation of an inner accretion disk ~a few tens of Rg. Time 
lags can also manifest in this series-of-discrete-interactions 
model (see Figure 5). Each stream self-interaction can produce 
the optical/UV emission similar to the shock-at-apocenter 
model. Furthermore, fluctuations from the last self-interaction 
—Just before the flow joins the inner accretion disk—can travel 
down to the black hole on a viscous timescale and modulate the 
X-rays. The viscous timescale in a thin accretion disk at a 
radius of a few tens of Rg is comparable to the observed lag 
(top left panel of Figure 4), and thus suggests that a TDF disk 
may be confined to within a few tens of Rg. 

We can also show that the expected properties of elliptical disks 
are consistent with the observed lags by constructing a simplistic 
model following G14. Recent numerical works have shown that 
the flow can remain elliptical for roughly 10 orbits (B17). Smaller 
mass black holes can maintain the ellipticity for longer duration. 
In such a model, we assume the emission follows the same radial 
profile as in a circular thin disk (Shakura & Sunyaev 1973), but 
the extent of the emitting region is limited by the size of the 
elliptical orbit of the most bound gas. Under such assumptions, 
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Figure 4. Top: the dependence of viscous timescale (top left) and the wavelength of emission (top right) on the radial distance in a Shakura & Sunyaev (1973) thin 
disk are shown. The observed lags are inconsistent with a circular thin disk. The solutions for different black hole mass, Eddington ratio, and a-viscosity parameters 
are shown. Bottom: viscous timescale vs. wavelength of emission for a standard thin disk. The observed time lags, shown as blue data points, are orders of magnitude 
faster than expected from a circular thin disk. A simplistic elliptical accretion disk model (G14) can reproduce the observed lags and is shown as a black curve. 


the radius dominating emission at a given frequency can be 
approximated as 


Rmax =l 
4r f. cm (max 1, a) Fan (v, Tae (R))AR 


Rmax n , 
anf x R(max[1. a) Fas(v, Re (R))dR 


where Rr is the tidal radius for 1 Mẹ star and Rmax is the 
apocenter radius for the most bound orbit. Fgg (v, Ter (R)) is the 
blackbody emission at frequency v and temperature Ter (R). 

To convert the radius at which emission takes place to the 
lags between given wavelengths and X-rays, we assume 
that the corresponding viscous timescale is given through 
tviscous = V- fg, where V is an ad hoc viscosity parameter and 
tee is the local free-fall time. In an elliptical disk, this ad hoc 
viscosity parameter is much higher than in a circular disk because 
we allow for larger “effective” a-viscosity and disk thickness. 

In Figure 4 (bottom panel), we plot the expected lags in the 
elliptical model as a function of wavelength obtained assuming 
Mpy, black hole mass = 1095 Mo, accretion rate = 0.2Mgqa, 
and Y = 0.05. The obtained relation fits the observed values of 


the lags very well. If the assumed value of Y = 0.05 is correct, 
the gas takes ~20 orbits to spiral down from given radius to the 
black hole, i.e., the gas would accreted quite rapidly toward the 
black hole. However, the dynamics of elliptical disks forming 
in tidal disruption events, and in particular their dissipation and 
emission profiles, are currently very poorly understood, and 
thus cannot be used to estimate the flow parameters 
unambiguously. 

In summary, we have presented a cross-correlation analysis 
that suggests that ASASSN-14li’s X-ray variations lag the 
optical/UV fluctuations by 32 + 4days. We propose that at 
least three models can explain these lags: (1) shocks-at- 
apocenter model, (2) an elliptical accretion disk model, and (3) 
series-of-discrete-interactions model. The shocks-at-apocenter 
model is similar to the model proposed by P15 but, in addition, 
we propose that the optical and the UV emission originates 
from physically distinct sites as opposed to the single 
photosphere scenario suggested by P15. ASASSN-14li’s data 
does not allow us to unambiguously detect the lags between the 
optical and the UV bands, but high-cadence optical and UV 
monitoring observations of a TDF would allow us to easily test 
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Figure 5. Schematic of the three likely stream self-interaction scenarios for the origin of the optical/UV emission from ASASSN-14li. In the shocks-at-apocenter 
model (top left), as clumps of debris trace the path shown, they first emit optical radiation, followed by UV and eventually—when near the black hole—modulate the 
X-rays. The time lag between the emission from these wavebands is dictated by the infall time between the optical, the UV, and the X-ray sites. In the elliptical disk 
model (top right), the time lags correspond to the local viscous timescale in the disk and are much faster compared to a circular disk of the same radial extent. In the 
series-of-successive-interactions model (bottom), fluctuations from the last self-interaction can propagate in the inner disk of a few tens of gravitational radii on the 


viscous timescale (see top left panel of Figure 4). 


this hypothesis and establish the origin of the optical/UV 
emission in other TDFs. 


This work is based on observations made with Swift, a 
mission that was managed and controlled by NASA's Goddard 
Space Flight Center (GSFC) in Greenbelt, MD, USA. AII the 
data used in the present article are publicly available through 
NASA's HEASARC archive. J.G. would like to thank Chelsea 
MacLeod, Anna Pancoast, and Yanfei Jiang for valuable 
discussions. Finally, D.R.P. would like to thank Poshak 
Gandhi, Simon White, and Margaret Trippe for valuable 
comments and suggestions. 
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